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ABSTRACT 

Narrow-band imaging of the rest-frame Lyman continuum (LyC) of galaxies at z ~ 3.1 has produced 
a large number of candidate LyC-emitting galaxies. These samples are contaminated by galaxies at 
lower redshift. To better understand LyC escape, we need an uncontaminated sample of galaxies that 
emit strongly in the LyC. Here we present deep Hubble imaging of five bright galaxies at z ~ 3.1 
that had previously been identified as candidate LyC-emitters with ground-based images. The WFC3 
F336W images probe the LyC of galaxies at z > 3.06 and provide an order-of-magnitude increase 
in spatial resolution over ground-based imaging. The non-ionizing UV images often show multiple 
galaxies (or components) within ~ 1" of the candidate LyC emission seen from the ground. In each 
case, only one of the components is emitting light in the F336W filter, which would indicate LyC escape 
if that component is at z > 3.06. We use Keck/NIRSPEC near-IR spectroscopy to measure redshifts 
of these components to distinguish LyC-emitters from foreground contamination. We find that two 
candidates are low redshift contaminants, one candidate had a previously misidentified redshift, and 
the other two cannot be confirmed as LyC-emitters. The level of contamination is consistent with 
previous estimates. For the galaxies with 2 > 3.06, we derive strong lcr limits on the relative escape 
fraction between 0.07 and 0.09. We still do not have a sample of definitive LyC-emitters, and a much 
larger study of low luminosity galaxies is required. The combination of high resolution imaging and 
deep spectroscopy is critical for distinguishing LyC-emitters from foreground contaminants. 

Subject headings: galaxies: high-redshift, galaxies: starburst, galaxies: intergalactic medium, ultravi¬ 
olet: galaxies 
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Surveys are now identifying large numbers of galax¬ 
ies at the reionization epoch (z ~ 7 — 8). Conse¬ 
quently, there has been much interest in how those galax¬ 
ies reionized the hydrogen in the intergalactic medium 
(IGM). Several studies have determined that the num¬ 
ber density of AGNs is not high enough to p r ovide 
the required ioni zi ng backgroun d (Ino ue e t all 120061 
Siana et all 120081: IWillott et all 120101 IMasters et al.l 
2012D . though some have suggested that faint active 
galactic nuclei (AGNs) may contribute significantly 
(|Glikman et al.l 120111: iFontanot et al.l I2012H . Because 
massive stars are the only other significant source of ion¬ 
izing photons, they are believed to be responsible for the 
reionization of intergalactic hydrogen. 

Whether or not stars reionized the intergalactic 
medium and provided the ionizing background for 1-2 
Gyr thereafter is contingent upon many uncertain pa¬ 
rameters: the total star formation rate density (SFRD) 
including the significant contribution from sources be¬ 
yond our current detection limits; the intrinsic ionizing 
spectrum of the galaxies; the fraction of the ionizing ra¬ 
diation produced by these galaxies that escapes into the 
intergalactic medium; and the “clumping” of the hydro¬ 
gen being ionized in the IGM at any epoch (and thus the 
recombination rate). 

Many efforts, both theoretical and observational, are 
being made to dete rmine the amount of undetected star 
formation (lAlavi et al.l 12 0 1 11) and the c lumping factor 
(|Pawlik et al.l l2009t iFinlator et all l2012f) . However, lit¬ 
tle is understood about the “escape fraction”, f esc , of 
ionizing radiation. Young stars form in dense molecular 
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clouds. In order for a significant fraction of ionizing radi¬ 
ation to escape into the IGM, lines of sight (LOSs) with 
low neutral hydrogen column-density must be cleared out 
or the stars must migrate out of the gas. Either process 
must occur within the short lifetime (< 10 Myr) of the 
O-stars that produce the ionizing radiation. 

It is not understood how the stars can be ex¬ 
posed on such short timescales, though several mecha¬ 
nisms have been proposed. These include supernovae 
winds and their effects on the Hi distribution in both 
the interstellar medium (ISM) and the circumgalactic 
medium (CGM) I Dove et all 120001 iClarke fc Oev 


iFuiita et al.1 120031 iRazoumov fc Sommer-Larsenl 1200 


intera ctions/mergers with other gala xies (IGnedin et al 


2008), and runaway massive stars dConrov fc Kratter 


20121) . Each of these proposed mechanisms has observa¬ 


tional signatures that would allow us to test which pro¬ 
cesses are responsible for high escape fractions and, thus, 
the reionization of the universe. 

Regardless of the precise mechanism, the first goal is to 
measure the average escape fraction and any dependen¬ 
cies it might have on galaxy properties or environment. 
Over the last 15 years, many attempts have been made 
to directly detect escaping Lyman continuum (LyC) 
radiation from galaxies of varying types and at various 


Steidel et al. 200l|; Giallongo et al. 2002|; 

Malkan et al. 

2003f Siana et al. 20071 

Crimes et al.l 2007 

; Cowie et al. 

20091: Grimes et al. 2009 

; Iwata et al.l 2009: Bridge et al. 


201ft ISiana et all 201 ()1 Vanzella et, al.l 2010fc : 

Nestor et al. 11201 lHLeitet et al.ll201lHBoutsia et al.ll201ll 


Vanzdla^_al. 20121 Nestor et al. 20131 IM ostard i et al.l 

20131: iLeitet et al.l 1201.31 iBorthakur et al.ll2014D . These 


studies must be conducted at 2 < 3.5 because the IGM 
becomes optically thick to ionizing photons at highe r 
redshift (jlnoue fc Iwatal 120081: iProchaska et al.1 [2009). 
Also, galaxies in the nearby universe must be receding 
at sufficient velocities such that a significant fraction 
of the ionizing spectrum is redshifted to wavelengths 
greater than_ 912 A to avoid absorption by Galactic 
H I dHurwitz et al.1 11997D . The most sensitive limits 
have been obtained at z ~ 1 with ultraviolet imaging 


or spectroscopy with 

the Hubble Space Telescope 

(HST) and GALEX 

Malkan et al.l 

2003 

Siana et al. 

2007 

: Cowie et al. 2009 

Bridge et al. 

2010; 

Siana et al. 

2010 

) and at z ~ 3 with optical data from the ground 

(Shaolev et al. 20061: Iwata et al.l 2009t Nestor et al. 

2011, 2013 Mostardi et a.1.1 2013^ The re 
studies at z ~ 1 and z ~ 3 are quite differer 
there are no detections from a sample of 

suits of the 
it. At z ~ 1, 
68 galaxies 
and a null 

(Siana et al. 2007, 2010; Bridge et al. 

2010 


At 2 : ~ 3 however, roughly 10% of galaxies h ave been 
detected with large ion i zing emissivities (|Iwata et al.1 
120091 iNestor et al.l l201lL I2013D . Because the limits at 
z ~ 1 are typically far deeper than the surveys at 2 ~ 3, 
the large number of detections at z ~ 3 suggests an 
averag e escape fraction that is significantly higher at 
z ~ 3 (ISiana et al.ll2010D . 

There are several possible explanations for the seem¬ 
ingly discrepant results at z ~ 1 and 2 ~ 3. First, the two 
samples were selected in different ways. In particular, 
most of the galaxies in the z ~ 3 samples were selected 


as Lya-emitters. Because both the Lya and LyC escape 
are dependent upon the distribution of neutral Hydro¬ 
gen, it is certainly possible that a Lya sel ectio n would 
also tend to find LyC emitters. However, INestor et all 
( 2011 ) saw lower Lya fluxes for galaxies with LyC de¬ 
tections, suggesting that Lya emission is not necessar¬ 
ily positively correlated with high LyC escape fractions. 
Second, the role of environment may be important. The 
large samples at z ~ 3 have been selected within large 
over densities to increase observing efficiency, but the in¬ 
creased star formation in these dense environments might 
affect the local ionizing background and the neutral hy¬ 
drogen columns in the circumgalactic medium. Finally, 
the z ~ 1 and z ~ 3 samples may be in different phases 
of their evolution. T he samples have similar SFRs and 
stellar masses (ISiana et al.l 12010 ). However, because the 
z ~ 1 samples were selected to be very luminous, they 
lie well above the star -forming main sequence at z ~ 1.3 
dWhitaker et al.l[2012D , and are likely in a transient star- 
burst phase. It is possible that the beginning stages of 
such a starburst may not be conducive to high LyC es¬ 
cape fractions (for example, because gas fractions are 
particularly high or because many are merging systems). 
Finally, it may be possible that other galaxy properties 
that have not been measured well ■ their sizes or circum¬ 
galactic medium, for example - can influence the LyC 
escape fraction. 

Many uncertainties remain in understanding ionizing 
photon escape and the differences in the 2 ~ 1 and z ~ 3 
studies. Foremost among these is the possibility of fore¬ 
ground contamination of the z ~ 3 LyC detected galax¬ 
ies. In some cases, it is possible that a faint star-forming 
galaxy lies along the line of sight to the target galaxy, 
but at significantly lower redshift. In this case a candi¬ 
date detection in the rest-frame LyC of the background 
galaxy could simply be non-ionizing emission from the 
foreground galaxy. At z ~ 1, the line of sight is much 
smaller and the flux densities required to probe signif¬ 
icant escape fractions are not as low. Therefore, there 
are far fewer foreground objects along the line of sight. 
In addition, the z ~ 1 studies have been done with the 
Hubble Space Telescope and its high spatial resolution 
allows for easier identification of foreground objects at 
small angular separation (or impact parameter). 

Some studies have used deep U band (which probes 
below the Lyman limit above 2 > 3) number counts 
to determine t he probability of f oreground contamina¬ 
tion at z ~ 3 (ISiana et ahl [20071: [Vanzella et al.l 1201 Obt 
INestor et al.112011k IMostardi et al.l I2013D and find that 
a significant fraction of LyC detections are in fact at¬ 
tributed to foreground contamination. This contamina¬ 
tion fraction is a strong function of the depth of the LyC 
images (since the number density of foreground objects 
rises steeply at fainter magnitudes, e.g. lAlavi et ahl 
l20li) . as well as the FWHM of the point-spread func¬ 
tion (as the PSF widens, contaminants at larger impact 
parameters cannot be recognized). 

These statistical corrections for foreground contamina¬ 
tion are useful when determining what fraction of galax¬ 
ies have detectable LyC escape, or for constraining the 
total ionizing emissivity of galaxies. However, it is useful 
to know for certain which galaxies truly have high escape 
fractions, as they need to be studied in greater detail to 
understand the mechanisms for escape. 
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There are two ways to mitigate the foreground contam¬ 
ination problem. First, deep spectroscopy of candidate 
LyC emitters may provide evidence of the foreground 
object if emission lines are detected at other redshifts. 
However, the foreground object may lack strong emission 
lines or the lines may be shifted out of the wavelength 
range covered by the spectrum. Other evidence for a 
foreground galaxy may be seen in absorption. Any fore¬ 
ground contaminant would lie at a very sm all impact pa¬ 
rameter (< 10 kpc), and lSteidel et al.! (1201 Ol ) have shown 
that galaxies generally exhibit strong metal-line absorp¬ 
tion at these distances. A high S/N detection of the 
continuum of the background galaxy is required to iden¬ 
tify foreground absorption and is very difficult to obtain 
for galaxies fainter than R = 24.5. 

Another way to minimize the frequency of foreground 
contamination is to increase the spatial resolution of the 
LyC imaging, as the contamination rate is proportional 
to the area of the seeing disk. If one resolves the back¬ 
ground galaxies (with half-light radii of ~ 0"2) instead 
of spreading the light over the typical seeing disc size 
(0"6—l'/O), the contamination rate goes down by a factor 
of ~ 9—25 and becomes almost negligible ([Vanzella et al.1 

120H . 


In addition to simply removing the possibility of fore¬ 
ground contamination, higher resolution images may al¬ 
low us to resolve the LyC-emitting regions of galaxies. If 
SNe winds clear low column density sightlines, the LyC 
may escape primarily from the regions with highest star 
formation surface density. However, if galaxy interac¬ 
tions or runaway stars are the primary mechanism, the 
LyC emission will be less concentrated than the non¬ 
ionizing UV continuum. 

We have set out to determine definitively whether can¬ 
didate LyC emitters have foreground contamination or 
indeed have high ionizing emissivities and, if the latter, 
which mechanisms allow high escape fractions. To do this 
we have followed up some of the brightest LyC-emitting 
candidates s e lected from ground-based observations by 
I Nest or et all (120111 ) in SSA22, a field with a large over¬ 
density of galaxies at z = 3.09 dSteidel et alJll998h . Be¬ 
cause of the large number of galaxies at the same red- 
shift, this field has been used to efficiently find candidate 
LyC emitters with a narrow-band filter just be l ow the ob¬ 
served wavelength of the Lyman limit ( Iwata et alJl2009t 
INestor et all [mill [mill . 


First, we obtained deep HST imaging of the rest-frame 
LyC of these galaxies to determine from which galaxy 
(or regions of the galaxy) the candidate LyC is escaping. 
Second, we obtained deep Keck near-IR spectroscopy to 
detect the rest-frame optical emission lines (primarily [O 
ill] A5007) at high spatial resolution (~ 0'/5) from the re¬ 
gions emitting apparent LyC emission to verify the red- 
shifts of the sources. The Hubble imaging allows us to 
identify possible contaminating faint foreground galax¬ 
ies, and the Keck spectroscopy allows us to identify the 
redshifts of all candidate LyC emitters. Together, these 
observations allow us to detect foreground contaminants 
and, therefore, identify the true LyC emitters. 

In Section [2l we discuss the Hubble observations. In 
Section [31 we discuss the Keck near-IR spectroscopic ob¬ 
servations. In Section [JJ we discuss each of the galaxies 
individually. Finally, in Section [5[ we put the findings in 
context with other recent work. 


2. HST OBSERVATIONS AND DATA REDUCTION 
2.1. Target Selection 

Our target galaxies are all bright, L > 0.5 L* LBGs 
at z > 3.06 in the S SA22 field and were originally cata¬ 
logued in iSteidel et al.1 (120031 ). The LyC region of these 
galaxies at z > 3.06 was imaged with a blue narrow- 
band filter with the Low Resolution Imaging Spectro¬ 
graph (LRIS) on Keck I (NB 3640, see Fig u re CD f or the 
narrow-band filter curves) by INestor et alJ (1201 111 . The 
galaxies were selected by their Lya emission at z ~ 3.09 
with a narrow-band filter (NB4980) or viaa Lyma n brea k 
selection using broad-band filters (ISteidel et, all 120031) . 
The Hubble WFC3/UVIS field-of-view (2.7' x 2.7') is con¬ 
siderably smaller than that of Keck/LRIS (5.5' x 7.6') 
and a very deep exposure is required to detect the rest- 
frame LyC of these galaxies with Hubble. Therefore, we 
chose to obtain a single, deep HST LyC image near the 
NW corner of the LRIS image that has several bright 
candidate LyC emitters. The bright galaxies lying in 
the single WFC3/UVIS pointing are MD32, C49, MD46 
and D17. There is an additional candidate, aug96M16, 
that does not lie in this Hubble field, but we obtained 
Keck/NIRSPEC spectroscopy of that source as well (see 
Section [3j. 


2.2. Hubble Observations 

We were allocated a total of 39 orbits in Cycle-17 (PID: 
11636). All primary observations are summarized in Ta¬ 
ble [D Most (32) of the orbits were dedicated to deep 
imaging with the F336W filter, which probes the rest- 
frame LyC of galaxies at z > 3.06 (see filter curve in Fig¬ 
ured). The required exposure time in the F336W filter is 
large, as the expected ionizing continuum is much lower 
than the non-ionizing UV continuum. Below ~ 4000 A, 
read noise is the dominant source of noise in UVIS imag¬ 
ing. Therefore, the F336W exposure times were long 
(half-orbit in length, 1325s) to minimize the number of 
read outs. The 32 orbits were split into two-orbit visits. 
In each visit, a standard four point dither (WFC3-UVIS- 
DITHER-BOX) was used to achieve sub-pixel sampling 
of the point-spread function (PSF). Each visit performed 
the same dither pattern, but at slightly different central 
pointings spread over ~ 3". 

In addition to the F336W images, we also obtained 
optical imaging (F606W, F814W, which probe the non¬ 
ionizing rest-frame UV continuum of the targets) with 
the Wide Field Camera on the Advanced Camera for Sur¬ 
veys (ACS/WFC) and near-infrared imaging (F110W, 
F160W, which probe the rest-frame optical contin¬ 
uum of the targets) with the WFC3 infrared channel 
(WFC3/IR). The field was imaged for one orbit in each 
of these additional filters with a standard four-point sub¬ 
pixel dither pattern (ACS-WFC-DITHER-BOX, WFC3- 
IR-DITHER-LINE). Two pointings were required with 
WFC3/IR to image all of the targets in the field because 
the IR channel has a smaller field-of-view than that of 
the UV channel. 

In addition, there is a 3-orbit ACS/WFC F814W image 
from a Cycle-13 program (PID 10405, PI: Scott Chap¬ 
man) which covers most of our field and which we incor¬ 
porated into our data. 
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Fig. 1.— The shaded region shows the composite spectrum of IShaplev et all (120031 ) of z ~ 3 LBGs, shifted to the redshift of the over 
density in SSA22, z = 3.09. The spectrum below rest-frame 920 A is just an extrapolation assuming a constant value in f u and multiplying 
by the average transmission of the intergalactic medium at 2 : = 3.09, expAlso plotted are the Keck/LRIS transmission curves 
for the NB4980 filter (used to find Lya-emitters at the redshift of the over density) and the NB3640 filter (used to probe the LyC photons 
from galaxies at 2 > 3.09). Finally, we plot the total system throughput for the WFC3/UVIS F336W filter - also used to detect LyC 
photons from galaxies at 2 : > 3.09. F336W is a broader and bluer filter than NB3640. Depending upon the line-of-sight IGM opacity, the 
F336W filter can detect more LyC photons than the NB3640. However, the average flux density across the filter will be much lower. 


TABLE 1 

HST Observation Summary 


Camera 

Filter 

Orbits 

Exptime 

Depth a 

WFC3/UVIS 

F336W 

32 

84800 

29.27 

ACS/WFC 

F606W 

1 

2140 

27.69 

ACS/WFC b 

F814W 

1 

2280 

27.28 

ACS/WFC b 

F814W 

3 

6144 

28.05 

WFC3/IR (East) 

F110W 

1 

2616 

28.07 

WFC3/IR (West) 

F110W 

1 

2616 

27.74 

WFC3/IR (East) 

F160W 

1 

2616 

27.42 

WFC3/IR (West) 

F160W 

1 

2616 

27.19 

WFC3/IR (aug96M16) c 

F160W 

1 

2612 

27.43 


a 3cr depth in a 0.5 // diameter aperture 

b Most of the area was observed previously with ACS/WFC F814W 
to three orbit depth in Cycle 13 (PID: 10405; PI: S. Chapman). 
c F160W image obtained in Cycle 17 (PID: 11735; PI: F. Man- 
nucci). 


Finally, one of the objects in this study (aug96M16) 
lies far outside the UVIS field-of-view. There is existing 
ACS F814W imaging (3 orbit depth, PID 10405) and 
F160W imaging (one orbit depth, PID 11735) that we 
reduced and shifted to the same astrometric system as 
the Keck/LRIS images. 


2.3. Data Reduction and Analysis 

The HST data were reduced using multidrizzle, the 
standard PyRAF routines provided by the Space Tele¬ 
scope Science Institute (STScI). The flattened frames 
were cleaned of cosmic rays and stacked onto images with 
matching WCS and aligned within a fraction of a pixel. 
Because sub-pixel dither patterns were used, we sample 
the PSF on a finer scale than the native pixel scales of 
0"04 for WFC3/UVIS, 0"05 for ACS/WFC and 0713 for 
WFC3/IR. Therefore we drizzled the data onto images 
with pixel sizes of 0"02, 0"04, and 0"08 for WFC/UVIS, 
ACS/WFC, WFC3/IR, respectively. Using stars in the 
field, we measure PSF FWHMs of 0707, 0709, and 0"18 
in the WFC3/UVIS, ACS/WFC, and WFC/IR channels. 


The output variance maps from multidrizzle were used 
to determine the depths of the images. These depths 
agree well with the pixel-to-pixel variations in the im¬ 
ages once t hey are corrected for co rrelated noise via the 
methods in iCasertano et all (I2000D . The 3cr depths are 
listed in Tabled] 

We combined the F110W and F160W images (weight¬ 
ing by the inverse variance in each image pixel) to pro¬ 
duce a deep detection imag e. We ran SourceExtrac- 
tor dBertin fe Arnoutsl 1 1 9961) to identify sources above 
the background fluctuations. The isophotes from the 
deep detection image were then mapped onto the finer 
pixel scales of the WFC3/UVIS and ACS/WFC im¬ 
ages. Because the images were drizzled to pixel scales of 
0702, 0704, and 0"08 (in WFC3/UVIS, ACS/WFC, and 
WFC3/IR, respectively), each pixel in the WFC3/IR im¬ 
age has either four ACS/WFC pixels or 16 WFC3/UVIS 
pixels in it. Therefore, the isophotes from the detection 
image can be mapped directly to the finer pixel images. 
Finally, we used our own code to determine isophotal 
fluxes in the same isophotes in all filters. The images of 
the targets are displayed in Figure [5] and fluxes are listed 
in Table [2] 

In addition to the fluxes, the estimated relative escape 
fractions (or the la limits) are listed as well. The rela¬ 
tive escape fraction, f^ c C re i is defined as in iSiana et al.l 
(120071) : 


f LyC 
J esc,rel 


(/l50o/ fLyc) stel 
(/l50o/ fLyC^obs 


exp (T H i t iGM) 


( 1 ) 


where /1500 and /l v c are the flux densities, in f v , at 
1500 A (F814W fluxes) and in the LyC (F336W fluxes, 
~ 820 A for these observations), respectively. The LyC 
flux density is assumed to be flat (constant) in f v be¬ 
fore IGM attenuation and the intrinsic amplitude of 
the Lyman break in th e ste llar SED is assumed to be 
(/ 1500 / fL y c) stei = 6.0 (ISiana et a l.l 120071) . We note that 
this ratio can be affected by short term variations in the 
star formation history. Specifically, after a rise in the star 
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formation rate, the LyC will be significantly enhanced 
compared to the 1500 A continuum. Conversely, shortly 
after a downturn in the star formation rate, the LyC 
flux can quickly decrease whereas the 1500Aflux will be 
largely unaffected. Using numerical simulations of galax¬ 
ies with realistic feedback prescriptions that broadly re¬ 
produce the sp read in the star-form i ng main sequence 
at the se epochs dWhitaker et al .1120121) , iDominguez et al.1 
(|2014f ) show that galaxies with M* > 10 9 Mq have a 
~ 0.09 dex (or ~ 23%) spread in the (/ 1500 / fL V c)stei 
ratio. 

In addition, the (/ 1500 / fL V c)stei ratio may be signif¬ 
icantly lower (by a factor of ~ 2) when including the 
effects of the binarity and rotation of massive stars, as 
the number of WR stars is increased and the main se- 
q uence lifeti me and effecti ve temperature are increased 
(jEldridee fe Stanwavll2009i ). 

The IGM transmission, exp(— tigm,hi), is integrated 
along the entire F336W filter bandpass. 


exp( tjgm,hi ) 


/ exp(-r A ) T\ dX 
f T\ dX 


( 2 ) 


where the IGM opacity, r, and the filter transmission, 
T, are both functio ns of wavelen gth. Using a similar 
technique as in lNestor et al.1 (12013H to simulate the IGM 
lines of sight, the value for the average IGM transmis¬ 
sion through the F336W filter curve at z = 3.09 was 
estimated to be 0.184 though the value can vary wildly 
among sight lines (16th and 84th percentile values are 
0.011 and 0.365, respectively). The most opaque sight 
lines can be ruled out if the NB3640 detections are real. 
The IGM transmission varies slowly with redshift. In 
Table m we list the average transmission value at each 
redshift and use that value to compute f^ c G re i- 


3. NIRSPEC OBSERVATIONS 

Near-IR spectra were obtaine d for four of the targ ets 
using the NIRSPEC instrument dMcLean et al.llT~998f ) on 
the Keck II telescope. We observed with NIRSPEC on 
18 and 19 August 2011. Conditions were photometric, 
with 0U5 seeing in the K band. The exposure times were 
6x900 seconds for C49 and MD32, 8x900 seconds for 
MD46, and 4x900 seconds for aug96M16. All targets 
were observed with a 0'.'76 x 42" long slit. As explained 
in Section 01 many of the targets are either elongated 
or have several detected clumps (or separate galaxies) 
within a 1"2 radius. The position angle of the NIRSPEC 
slit (listed in Table 0]) was chosen in the direction of elon¬ 
gation or to collect spectra of multiple objects. Since 
the primary spectral features of interest were [OIII] A5007 
and H/3 at z ~ 3.06 — 3.30, we used the N6 filter, a broad 
H + K filter centered at 1.925/mi with a bandwidth of 
0.75/im. The spectral resolution of our NIRSPEC ob¬ 
servations was ~ 15 A, as determined from sky lines 
measured in the N6 filter. This resolution corresponds 
to Av ~230 km s -1 (R~1300). Due to the faint nature 
of these objects in the K band, we acquired each target 
using blind offsets from a bright star in the surrounding 
field. We returned to the offset star between each inte¬ 
gration of the science target to recenter and dither along 


the slit. In all cases we dithered back and forth between 
two positions near the center of the slit. 

3.1. NIRSPEC Data Reduction 

Data reduction was performed using a method simi¬ 
lar to that described in iLiu et all (1200811 , where the sky 
background was subtracted from the two-dimensional un¬ 
rectified science images using an optimal method dKelsonl 
120031 G. D. Becker 2006, private communication). The 
sky lines were fit with a low-order polynomial and a b- 
spline fit was used in the dispersion direction. After back¬ 
ground subtraction, cosmic rays and bad pixels were re¬ 
moved from each exposure. The individual images were 
then rotated, cut out along the slit, and rectified in two 
dimensions to take out the curvature both in the wave¬ 
length and spatial directions. The final rectified two- 
dimensional exposures for each object were then regis¬ 
tered and combined into one spectrum. 

Each of our targets had previously been spectroscopi¬ 
cally confirmed as a z ~ 3 galaxy based on optical (i.e., 
rest-frame UV) spectroscopy. Our goal was to determine 
if the specific regions associated with apparent LyC emis¬ 
sion were also at z > 3.06, high enough redshift that the 
NB3640 and F336W filters provided a clean probe of the 
LyC spectral region. Accordingly, in addition to confirm¬ 
ing the redshift for the main target, we also searched the 
NIRSPEC slit for additional emission-line features offset 
from the main region of non-ionizing UV flux, and per¬ 
haps originating from the same regions as the NB3640 or 
F336W flux. 

In the spectra of C49, aug96M16, and MD32, we dis¬ 
covered evidence for an additional emission-line spec¬ 
trum coincident with or very near the location of 
NB3640/F336W emission (referred to, respectively, as 
C49JNTW, aug96M16_E & aug96M16_W, and MD32.NE), 
while MD46 yielded only a single detected spectrum cen¬ 
tered on the location of the non-ionizing UV continuum. 

One-dimensional spectra were extracted from the two- 
dimensional reduced image at the location of any de¬ 
tected emission lines. The corresponding lcr error spec¬ 
tra were also extracted. The average aperture size along 
the slit was l."3, with a range of l."l to 2."0. Each 
one-dimensional spectrum was then flux-calibrated us¬ 
ing A-type stars acc ordi ng to the metho d described in 
iShaplev et al.l (1200511 and lErb et al.l (12003H and placed in 
a vacuum, heliocentric frame. The final flux-calibrated, 
one-dimensional spectra are plotted in Figure 0] 

3.2. NIRSPEC Measurements 

In the rest-frame optical, the strongest features are 
[OIIIJA5007 and H/3. We determined the [OIIIJA5007 
centroid, flux and FWHM, by fitting a Gaussian profile 
to each emission line using the IRAF task, splot. The 
nebular velocity dispersion was calculated as oy(0[III]) 
= (FWHM/2.355) x(c/A 0 t, s ), where A 0 bs is the observed 
wavelength of [OIIIJA5007. The FWHM was determined 
from the subtraction of the instrumental FWHM (esti¬ 
mated from the widths of night sky lines) from the ob¬ 
served FWHM in quadrature. Uncertainties in the veloc¬ 
ity dispersion were det ermine d using the same method 
described in lErb et al.l (l2006bih Where detected, we also 
measured the properties of H/3 emission. 

A Monte Carlo approach was used to measure the 
uncertainties in the emission-line centroid, flux, and 
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NB3640 (LyC) F336W (LyC) F606W F814W F110W F160W 



Fig. 2.— Hubble imaging of the five bright LBGs with Keck NB3640 detections from I Nest or et ahl Q2011I) (candidate LyC-emitters). The 
stamps are 3" on a side, with tick marks at 0 / . / 5 intervals (the aug96ml6 stamp at the bottom is 3" on a side). In C49, MD46, and MD32, 
the F336W (possible LyC) flux is emitted from only one of multiple clumps. The D17 detection is the low surface brightness detection in 
the top left of the NB3640 image. The HST imaging is not as sensitive to low surface brightness emission so it can not definitively detect 
it. In the aug96ml6 stamp the central galaxy has strong emission in the NB3640 filter but this galaxy lies outside of the footprint of the 
UVIS pointing. 


TABLE 2 

Galaxy Properties 


Galaxy z spec exp (~t igm Y f^° rel bc F336W bd F606W d F814W d F110W d F160W d 


C49 

3.153 

0.158 

< 

0.088 


< 

0.037 

5.74 

± 

0.16 

6.37 

± 

0.10 

6.69 

± 

0.14 

8.13 

± 

0.24 

C49.NW 

2.029 




0.40 

± 

0.026 

1.73 

± 

0.19 

1.28 

± 

0.07 

1.89 

± 

0.10 

3.28 

± 

0.17 

MD46 

3.091 

0.0183 

< 

0.066 


< 

0.060 

9.31 

± 

0.25 

11.9 

± 

0.16 

11.1 

± 

0.18 

14.4 

± 

0.33 

MD46.NW 




3.62 e 

0.40 

± 

0.022 

0.87 

± 

0.09 

1.45 

± 

0.06 

1.47 

± 

0.10 

1.23 

± 

0.18 

MD32 

3.102 

0.178 

< 

0.451 


< 

0.025 

0.41 

± 

0.11 

0.75 

± 

0.07 

0.78 

± 

0.07 

1.43 

± 

0.12 

MD32.NE 

2.885? 1.964? 




0.56 

± 

0.038 

1.74 

± 

0.16 

2.22 

± 

0.10 

3.28 

± 

0.11 

5.41 

± 

0.20 

aug96M16_E 

3.095 










2.73 

± 

0.17 




6.31 

± 

0.45 

aug96M16 











3.87 

± 

0.24 




9.04 

± 

0.49 

aug96M16_W 

3.291 










3.22 

± 

0.14 




4.06 

± 

0.32 

D17 

3.069 

0.194 

< 

0.081 


< 

0.031 

4.10 

± 

0.14 

4.78 

± 

0.20 

4.61 

± 

0.13 

6.37 

± 

0.22 


a Defined in Equation [2] Here we use the average transmission of 500 simulated sightlines. 
b “<” denotes la limits. 
c Defined in Equation [I] 

d Flux density in units of 10 —30 erg s —1 cm -2 Hz -1 . 

e For this calculation, MD46.NW was assumed to have the same redshift as MD46, z = 3.091, and the same value of transmission 
through the IGM. An f esc re i > 1 means that the measured flux is larger than assumed. This could be caused by a very recent 
burst in star formation ( Dominguez et aLll2014j) . The number may be artificially high if line of sight transmission is higher than 
average. 
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C49 z = 3.1532 



C49_NW z = 2.0286 




bfl MD32 z = 3.0968 


S-. 

CD 




40 
20 
0 

-20 

1.98 2 2.02 2.04 2.06 2.0B 



40 
20 
0 

-20 

1.96 1.98 2 2.02 2.04 2.06 


MD46_NW 

-|—i—i—i—|—i—i—i—|—i—i—i—|—i—i—i—|—i—i—i—|—i— r~ 


l.-l—l.-t-l-L -« -4-.t-.l- A.. 


MD32_NE z = 2.8844 



aug96ml6_W z = 3.2899 


aug96ml6_E z = 3.0938 




Wavelength (/xm) 


Fig. 3.— The ID NIRSPEC spectra of all components from which line emission was detected. The lcr error spectrum is plotted as 
a dotted curve, and offset vertically by 20 x 10“ 19 ergs s —1 cm -2 A -1 from the corresponding science spectrum in each panel. The 
associated redshifts are given above each spectrum. The candidate LyC emitter, C49_NW, is found to be foreground emission from a 
galaxy at z = 2.0286. The candidate LyC emitter, MD46_NW has no detectable emission at the redshift of MD46 (z ~ 3.09, or any other 
redshift) so it cannot be verified as a LyC emitter based on these data alone. The candidate LyC emitter, MD32_NE, is found to have 
a redshift of z = 2.8844 if the detected line is [O ill] A5007. For aug96ml6, the NIRSPEC spectrum show us that our initial redshift was 
assigned to the incorrect galaxy. The central galaxy, which exhibits emission in the Keck NB3640 filter is not at z=3.2899. No redshift has 
been definitively identified for the candidate LyC emitter. 
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FWHM. For each object 500 simulated spectra were cre¬ 
ated by perturbing the flux at each wavelength of the 
true spectrum by a Gaussian random number with the 
standard deviation set by the level of the the la error 
spectrum. Line measurements were obtained from the 
fake spectra in the same manner as the actual data. The 
standard deviation of the distribution of measurements 
from the artificial spectra was adopted as the error on 
each centroid, flux, and FWHM value. Upper limits on 
the H/3 line fluxes were obtained for non-detections by 
assuming the same FWHM as the [OIII]A5007 emission 
line and performing the same Monte Carlo measurements 
as described above. Based on [OIII]A5007 emission line 
centroid measurements, we obtained the systemic red- 
shifts for each object. In one case (MD32_NE) the single 
detected emission line was assumed to be [O III] A5007. 
If this line identification is correct, [Oiii]A4959 and H/3 
would not be covered by wavelength range of that spec¬ 
trum. All NIRSPEC measurements are listed in Tabic [3] 

4. RESULTS 

Because there are only five galaxies in our sample and 
each has different properties, we have chosen to discuss 
the objects individually. We first discuss the HST imag¬ 
ing, as the resolved images informed the strategy for the 
near-IR spectroscopy (slit orientation and required see¬ 
ing for each galaxy). In all cases one must keep in mind 
that the Keck narrow-band filter (NB3640) used to iden¬ 
tify the LyC emission has a width corresponding to ~ 25 
physical Mpc, significantly smaller than the mean free 
path of a photon at the Lyman limit at z ~ 3.09 (~ 100 
Mnc. IProchaska et all 120091 : lO’Meara et al.ll2012l l. The 
HSTF336W filter, which we have used to detect the LyC 
of these galaxies, has a similar red cutoff (near the Lyman 
limit) but is a factor of ~ 5.5 times wider, with the addi¬ 
tional transmission blueward of the NB3640 filter (Figure 
[T]). At the bluer wavelengths probed with the F336W fil¬ 
ter, but not probed by the NB3640 filter, the average 
line-of-sight through the IGM toward a z = 3.09 galaxy 
will have significantly less transmission. Thus, the typi¬ 
cal galaxy at z = 3.09 should appear significantly fainter 
(1.0 magnitudes) in the F336W filter than in the Keck 
NB3640 filter (see Section [5] for more details). 

4.1. C49 

The galaxy C49 was first identified as a LyC emit¬ 
ter in deep rest -frame ultraviolet spectroscopy from 
Keck fehaplcv et al.ll200fih . and again with narrow-band 
imaging just below the Lyman limit (llwata et alJl2009t 
INestor et al.l 1201 0) . At ground-based reso lution (0"8), 
there was already a suggestion in Figure 3 of lNestor et al.1 
(120111 1 that the LyC flux was only emitted from the NW 
corner of the non-ionizing UV-emitting area of C49. 

The HST imaging resolves C49 into two sources sepa¬ 
rated by 0'.'65 (5.0 proper kpc at z = 3.153) in all filters. 
The fact that the two clumps are distinct at rest-frame 
optical wavelengths suggests that these are two separate 
galaxies rather than two UV-bright star-forming regions 
of the same galaxy. Because the galaxy to the SE is 
brighter and has a previously known redshift (z = 3.153), 
we are referring to it with th e original C49 label and 
the nearby galaxy as C49_NW. IShaplev et al.l (|2006l) had 
previously obtained a deep optical spectrum (rest-frame 
UV) of this system, but the two sources were unresolved 


and most of the light from both objects fell within the 
1"2 wide slit. Because C49 is five times brighter than 
C49_NW in the optical, it is the galaxy with the previ¬ 
ously reported redshift. 

The F336W (LyC) image shows strong emission from 
C49_NW, with no detectable emission from C49. We 
derive a la limit of f eS c,rei < 0.088 (lcr) for C49. If 
C49_NW is at the same redshift as C49, or at any other 
redshift above z = 3.06, then the F336W flux would be 
a direct detection of LyC emission. 

The NIRSPEC slit was oriented to lie along both ob¬ 
jects (PA = —38°). The seeing was sufficient (0"5) to 
marginally resolve [O ill] lines from both components if 
they were both at the same redshift. The 2D spectrum 
is shown in Figure [4] The strong [Om] AA4959, 5007 
emission lines are clearly detected at the redshift of C49. 
There is no indication of fainter [Om] emission at the 
same redshift from the companion. Instead, an emission 
line is detected at shorter wavelength (A = 1.988 /an) 
- 0"5 from C49 in the direction of C49.NW. The two 
most likely identification of this line are [Oiii]A5007 or 
Ha, indicating redshifts of z = 2.9698 or z = 2.0286, 
respectively. 

Given that C49_NW is likely at one of these two red- 
shi fts, we re-ex amined the deep LRIS spectrum of C49 
(IShaplev et al.l[2006h to look for absorption features from 
a foreground source. As mentioned above, the impact 
parameter is ~ 5 kpc, and strong absorption from both 
the interstellar and circumg alact ic medium of C49_NW 
would be expected USteidel et al.l 1 2 01(1) . In Figure [5] we 
plot the deep Keck/LRIS spectrum. There are several 
strong absorption features associated with z = 2.0286 in 
the Lya forest (Si ii A1260, O i / Si n AA1302,1304, C ii 
A 1334, Si ii A1526) and at the red end of the spectrum (Fe 
ii AAA2343, 2374, 2382). Therefore, we conclude that the 
emission line in the NIRSPEC spectrum is Ha a nd the 
flux detected at A 5. 378 0 A in both lShanlev et ahl (I2006D 
and INestor et al.1 (|201 lh is not LyC emission. The fore¬ 
ground contamination was not originally recognized with 
the LRIS spectrum alone because all of the strong UV ab¬ 
sorption features are in the Lya forest of C49 (except for 
the faint Fe ii lines at the red end of the spectrum where 
the sensitivity is poor). The additional detection of an 
emission line in the rest-frame optical allowed us to iden¬ 
tify the absorption features as associated with C49_NW 
and not with foreground Lya lines. 

4.2. MD46 

MD46 is a bright LBG at the redshift of the overden- 
sity (z — 3.09) and wa s iden tified as a candidate LyC 
emitter by Iwata et al.l (|2009T ) and INestor et al.l (|2011jh 
llwata etalT l 20091 ) suggested that the LyC was offset 
from the brightest region of MD46. The HST images 
of the rest-frame non-ionizing UV and optical contin¬ 
uum show several clumps. As seen in Figure [2] most of 
the emission is coming from two bright clumps and two 
fainter clumps all within a O'!7 diameter region, which we 
are collectively calling MD46. It is this bright cluster of 
objects that has the known spectroscopic redshift. There 
is an additional object 1"2 (9.2 kpc) to the NW that we 
will refer to as MD46JNW. MD46 is not detected in the 
HST F336W (LyC) image but MD46_NW is strongly de- 
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TABLE 3 

NIRSPEC Measurements 


Object 

Slit PA a 

2 b 
^em 

Ha Flux c 

Ha FWHM d 

[OIII] Flux c 

[OIII] FWHM d 

M[OHI]) e 

H/3 Flux c 

C49 f 

-38 

3.1532 



7.5 ± 0.25 

17 ± 0.7 

54.7+« 

<1.4 

C49.NW 

-38 

2.0286 

1.4 ± 0.16 

14 ± 2.6 





MD46 

-34 

3.0857 



5.6 ± 0.15 

16 ± 0.5 

27.11]° 

1.4 ± 0.15 

MD46_NW f 

-34 




<0.40 



MD32 f 

37 

3.0968 



1.2 ± 0.14 

14 ± 1.4 


<0.36 

MD32.NE 

37 

2.8844*5 



2.3 ± 0.19 

15 ± 1.7 



aug96ml6_W f 

90 

3.2899 



3.0 ± 0.24 

15 ± 1.8 

<50.48 

<0.82 

aug96ml6_E 

90 

3.0938 



2.9 ± 0.34 

18 ± 2.1 

60.2l^ 

0.8 ± 0.25 

aug96m!6 f 

90 




<1.2 





a Slit position angle in degrees East of North. 

b Rest-frame optical nebular emission-line redshift measured from [OIIIJA5007 or Ho;. 
c Emission-line flux and error in units of 10“ 17 ergs s —1 cm -2 . 

d Observed FWHM and error in units of A. The FWHM is a raw value, uncorrected for instrumental broadening. 
e Nebular velocity dispersion measured in units of km s —1 . 

f Upper limits for [OIII] and H/3 fluxes are quoted at the 3 cr level. For aug96M16 and MD46.NW, limits were calculated by 
assuming the same redshift, respectively as aug96M16_E (z = 3.0938) and MD46 (z = 3.0857). 

s Redshift estimated assuming that the single emission line detected corresponds to redshifted [OIIIJA5007. If the emission 
line is instead redshifted Ha, the redshift MD32_NE is z = 1.9635. 






Fig. 4. The 2D NIRPSEC spectra (right) and the HST F814W stamps of the targets (left) rotated so that the NIRSPEC slit is oriented 
up-down, to match the 2D spectra. The y-direction spans 3" for the top three galaxies, 3" for the bottom galaxy and the spectra and 
images have the same spatial scale. The red dashed vertical line denotes the wavelength above which the [Om]A5007 line must lie in order 
for the NB3640 and F336W filters to measure purely LyC, without contamination from flux at wavelengths above the Lyman limit. Both 
C49 and MD32 show emission lines at two redshifts with the lower redshift lines coming from the higher (in the y direction) of the two 
components (which are the components detected in NB3640 and F336W). In MD46 there is no detection of emission lines from the faint 
companion which is detected in NB3640 and F336W. Even if the galaxy were at the same redshift, we would be able to resolve [Oin]A5007 
emission from it because the separation of the two components (l / . / 2) is much larger than the seeing (O'.'5). Finally, we detect emission 
lines for two of the galaxies in the region of aug96M16 but, given the vertical separation of the lines, they are associated with the upper 
and lower galaxies. We have not identified a redshift for the central galaxy, which is the one of interest as it is detected in the NB3640 and 
F336W filters. 


tected. We derive a la limit of f eS c,rei < 0.066 (ler) 
for MD46. If MD46_NW has the same redshift as the 
brighter counterpart to the SE (or any redshift above 
0 = 3.06), then the detected F336W flux is LyC emis¬ 
sion. 

As seen in Figure [2] the Lya emission is clearly emit¬ 
ted from the NW side of MD46, in the direction of 
MD46JNW, pos sibly indicati n g that MD46JMW is at the 
same redshift. iNestor et al.1 (1201311 use the Lya image 
and a newly-obtained deep rest-frame UV spectrum cov¬ 
ering Lya showing that the Lya emission extends from 
MD46 in the N and W directions and coincides with 
MD46_NW. They therefore conclude that MD46_NW is 
likely at the same redshift as the brighter MD46. Ad¬ 
ditional evidence for this interpretation is that there are 
no unexplained strong absorption lines in the high S/N 
rest-UV spectrum of MD46, which would arise from pas¬ 
sage through the halo (within ~ 10 kpc) of MD46_NW 


if it were in the foreground. However, neither of these 
arguments (the extended Lya emission or the lack of ab¬ 
sorption in the MD46 spectrum) definitively proves that 
MD46_NW is at the same redshift as MD46. First, as 
seen in C49, if the foreground object is at the right red¬ 
shift (z ~ 2), most strong absorption features will either 
be in the Lya forest or fall off the red end of the spec¬ 
trum. Also, because Lya is a resonance transition, Lya 
emission is commonly seen extended beyond the stellar 
extent (and often asymmetric). Thus, a definitive confir¬ 
mation must come from a non-resonance, emission line 
such as [O ill] A5007. 

The NIRSPEC slit was oriented to obtain spectra of 
both objects, MD46 and MD46_NW (PA = -34°). The 
separation between the two components (1"2), is large 
enough that emission lines at the same wavelength from 
the two targets would easily be resolved if at the same 
redshift. [O iii]AA4959,5007 and H (3 from MD46 are de- 
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Fig. 5.— LRIS spectrum of C49 (and C49_NW in the same slit) from [i^haplev et al.l <120061) . Red dashed lines denote the wavelengths of 
strong UV absorption lines at the redshift of C49 (z = 3.1532, strongly detected lines labelled in red at the bottom of the plot). The blue 
dashed lines denote those same lines at z — 2.0286 (strongly detected lines labelled in blue at the top of the plot). Within the Lya forest 
at A < 5049 A(denoted with a thick black dashed line), we still see strong absorption coinciding with Si II A1260, O I / Si II AA1302,1304, 
C II A1334, Si II A1526, and Fe II AAA2343, 2374, 2382. Therefore, we conclude that the emission line detected in the NIRSPEC spectrum 
for C49_NW is Ha, and C49_NW has z=2.0286. The redshift of C49_NW would have been difficult to determine with this LRIS spectrum 
alone because all of the strong absorption features are in the Lya forest or at the red (low S/N) edge of the spectrum. 


tected at high significance (S/N = 37 in the case of [O 
m]A5007). However, no emission lines are seen at the lo¬ 
cation of MD46_NW, with a 3er upper limit of 4.0 x 10~ 18 
erg s -1 cm -2 . 

We can estimate the expected [O ill] A5007 emission 
line flux for MD46_NW to determine whether our obser¬ 
vations are sensitive enough to detect the line if it lies 
at the same redshift as MD46. We make the assump¬ 
tion that the ratio of the [O ill] A5007 lines of the two 
galaxies is the same as the ratio of the rest-frame UV 
fluxes (as seen in the F814W filter). This is only true 
if the [O ill] emission traces the star formation rate and 
if there is similar excitation and dust extinction in both 
galaxies. The ratio of F814W fluxes is 8.21, giving a pre¬ 
dicted [O hi] A5007 flux for MD.NW of 6.8 x 10 -18 ergs 
s -1 cm -2 , which would be a 5.Id detection. Thus, the 
lack of detection suggests that MD46_NW might be at a 
different redshift such that its strong lines lie outside of 
the K-band (or on a sky line). Of course, this argument 
assumes that the [O ill] emission is proportional to the 
rest-UV continuum emission, but the T\om]/L jjv ratio 
is a function of star formation history, metallicity, and 
dust extinction. One might expect that the lower mass 
galaxy (MD46_NW) has a lower metallicity and, th us, a 
higher [O III] to Balme r line ratio (eg. lErb et al1 f2006a: 
iDomfnguez et al .1120131) . indicating that its [O III] to star 
formation ratio should be high. However, MD46JMW is 
redder in F606W-F814W color, suggesting that it may 
be dustier. Unfortunately, given the dispersion in UV- 
to-[0 ill] flux ratios seen in star-forming galaxies, it is im¬ 
possible to definitively determine with these data alone 
whether or not MD46_NW is at the same redshift as 
MD46 and, therefore, whether or not it is a LyC emitter. 

4.3. MD32 

MD32 is the faintest target in our sample and is an 
elongated galaxy with an extent of ~ 1"2, which ap¬ 
pears to be edge on. The F814W image shows two bright 


clumps to the NE and SW and a fainter, smaller clump 
between them (~ 0"2 of space between the two regions). 
The NE region is brighter and strongly detected in the 
F336W (LyC), as is the faint clump near the center. The 
SW clump is undetected in the F336W image. 

The NIRSPEC slit was placed in the direction of elon¬ 
gation of MD32 (PA = 37°) to search for emission from 
the two components or to identify rotation in what could 
be an edge-on disk. Faint line emission at the expected 
wavelength of [OIIIJA5007 at 2 = 3.0968 is detected, con¬ 
firming the redshift determined via identifi catio n of Lya 
emiss ion detected in the LRIS spectrum ([Nestor et al.1 
I20l3h . We also detect a stronger emission line at shorter 
wavelengths (A = 1.945 /nn). This line is offset nearly 1" 
along the slit in the NE direction, corresponding to the 
region of MD32 detected in the F336W image. There¬ 
fore, we conclude that the F336W-detected portion is at 
a lower redshift (z = 2.8844 or z = 1.9635 if the line 
is [O ill]A5007 or Ha, respectively), and is therefore not 
escaping LyC. Because MD32 has long been identified as 
being at z = 3.09, we refer to the SW clump as MD32 
and the NE clump as MD32_NE. We derive a Id limit of 
fesc,rei < 0.451 (lcr) for MD32. 

The importance of high-resolution imaging combined 
with rest-frame optical spectroscopy is demonstrated well 
here. The HST F336W image allowed us to determine 
definitively from which clump the NB3640 emission was 
emitted, and the rest-frame optical spectrum determined 
that the NB3640-detected clump is at lower redshift. 

4.4. aug96M16 

The object aug96M16 is perhaps the most complex of 
our targets as there are three galaxies lined up in the E- 
W direction. aug96M16 does not fall within the F336W 
footprint, but we have reduced existing HST F814W and 
F160W imaging of this system. The center and eastern 
galaxies are confused at ground-based resolution, but the 
HST images reveal that there are three distinct galaxies. 
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We will refer to these objects, from East to West, as 
aug96M16_E, aug96M16, and aug96M16_W. 

The centroid of the NB3640 coincides with the middle 
of the three galaxies, aug96M16. This gala xy was erro¬ 
neous ly as signed a redshift of z = 3.29 in INestor et al.l 
(|2011f) but INestor et al.l (|2013f ) corrected this and iden¬ 
tified the westernmost galaxy, aug96M16_W as having 
this spectroscopic redshift. Thus, aug96M16 should not 
have been included as a candidate LyC emitter because 
its redshift was unknown. Our NIRSPEC spectra were 
obtained before this error was realized. Nonetheless, if 
the central galaxy is at z > 3.06, the detected NB3640 
emission would be LyC emission. 

The NIRSPEC slit was aligned E-W (PA= 90°) so 
that all three objects fell on the slit. Two sets of 
[O m]A5007,A4959 emission lines are detected in the 
2D spectrum, corresponding to redshifts of 3.0938 and 
3.2899. T he latter reds hift ag rees with the redshift re¬ 
ported in INestor et al.l (120111 ). If the galaxy emitting 
flux in the NB3640 were at either of these redshifts, then 
the NB3640 would be LyC emission. Because we don’t 
know the absolute position of the slit on the sky, we use 
the distance between the objects to discern which lines 
correspond to which galaxies. The physical separation 
along the slit is 3" 24. This corresponds to the separa¬ 
tion between the easternmost and westernmost galaxy. 
Thus, the westernmost galaxy has a spectr o scopic red¬ 
shift of 3.2899 (consistent with INestor et al.l 12013 ) and 
the easternmost galaxy has a spectroscopic redshift of 
z=3.0938. We did not detect any strong emission lines 
from the central galaxy that is detected in the NB3640 
filter. Thus, the central galaxy can not be confirmed as 
a definitive LyC emitter. 

4.5. D17 

D17 is a bright compact galaxy located in the proto¬ 
cluster at z = 3.09, and was identified as h aving dif¬ 
fuse L yC emission offset ~ 1" to the NE by INestor et alJ 
{2CQj3). There is not a galaxy detected at the location 
of the emission at longer wavelengths. However, it is 
physically possible to have nebular LyC emission from 
a cloud that is photo-ioniz ed by LyC escaping from a 
nearby galaxy (jlnouel l2010h . However, such models re¬ 
quire extreme stellar p opulations (< 1/50 Zq, age < 1 
Mvr. llnoue et aLlfeOllh . and one of the best cases for such 
nebular LyC emission (object ‘a’ in llnoue et all 1201 ll) 
was found to actually be due to a foreground galaxy con¬ 
taminating t he spe ctral energy distribution (LAE 003 in 
INestor et al.ll2013j ). 

Because D17 lies in our F336W image footprint, we in¬ 
vestigated whether we also detect this diffuse emission in 
our image. The reported NB3640 emission near D17 was 
just above the detection limit, with an AB magnitude 
of 27.0. Assuming an average line-of-sight through the 
IGM and a flat (in /„) LyC spectrum, we would expect 
the F336W magnitude to be 28.0 (AB). At 1 arcsec 2 , 
that flux level corresponds to 4.3er. There is no signifi¬ 
cant detection in the F336W image. However, because 
the IGM transmission can vary significantly with wave¬ 
length (it can be transparent at 3640 A and opaque at 
shorter wavelengths), the F336W non-detection does not 
invalidate the detection in NB3640. Therefore, our image 
can not confirm the existence of offset nebular continuum 


NE of D17. 

We did however use our F336W image to estimate the 
LyC escape fraction “down-the-barrel” toward D17. Us¬ 
ing an aperture defined by the D17 near-IR image, we 
derive a la limit of f e sc,rei < 0.081. 

Because there are not other galaxies around D17 and 
the one galaxy has a known redshift, a NIRSPEC spec¬ 
trum was not obtained for this system. 

5. DISCUSSION 

There have been several surveys for LyC-emitters at 
this redshift, some detecting few or no LyC-emitters 
(|Boutsia et al.l [20TTI : lYanzella et al.l 120121). some claim¬ 
ing dozens of candidates (llwata et alJl2009t INestor et al.l 
1201 ll 120131 : IMostardi et al.ll2013f) . It can be confusing to 
synthesize these results into a unifying picture. As we 
have seen in this paper, it is critical to carefully consider 
contamination and its effects on each survey when trying 
to compare results. 

The rate of contamination from foreground galaxies is 
quite low at these wavelengths. As an example, the odds 
of a foreground galaxy with Uab < 27.5 lying within 
1"0 of a target galaxy at z ~ 3 is roughly 8.5% given 
the completene s s-corre cted surface densit y tabulated in 
lYanzella et al.l |2010b) and taken from iNonino et al.l 
(2009). But detection of LyC emitters is also rare. Ob¬ 
viously, if the detection rate is near the contamination 
rate, then most or all of the candidate LyC emitters will 
be spurious. It is therefore critical that any survey for 
LyC-emitting galaxies accurately calculates the contam¬ 
ination rate and has a significantly higher detection rate 
if it is to claim a bona-fide LyC detection. We list here 
several factors that can affect the contamination rates 
and detection rates in these surveys. First, the factors 
affecting the contamination rate, man y of which have 
been outlined before b y other authors (|Siana et alJl2007t 
lYanzella et al.|[2010bl) . 

• Depth: With increasing depth, fainte r and more 
nume rous foreground galaxies (e.g. lAlavi et al.l 
120141 can be detected. 

• Spatial resolution: At lower resolution, foreground 
galaxies at larger angular separation from the tar¬ 
get galaxy can contaminate the photometry. 

• Redshift: At higher redshift, because the line-of- 
sight to the target galaxy is much larger, the vol¬ 
ume of the contaminating cylinder is larger and the 
number of contaminants will therefore be larger. 
This is part of the reason that contamination con¬ 
cerns are much smaller for t he z ~ 1 samples of 
iSiana et~akl (|2007l . f20Toh and lBridee et al.l (|2010lh 

Generally, one can observe the surface density of ob¬ 
jects in the field at the wavelengths of the LyC obser¬ 
vations, and determine the rate at which an interloper 
would lie in the the seeing disk. The example above, 
because it is very deep and uses a fairly large, 1" ra¬ 
dius aperture, results in an 8.5% contamination rate. Of 
course, one also has to consider the flux distribution of 
expected interlopers as well. 

Below are the factors affecting the detection rate of 
galaxies. 
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• Depth: Deeper data will be sensitive to lower 
escape fractions, increasing the detection rate. 
Deeper data are also sensitive to escaping LyC from 
less luminous galaxies. 

• Redshift and IGM opacity: This factor is often 
overlooked. When looking at rest-frame wave¬ 
lengths near the Lyman limit (e.g. A res t = 850 — 
910 A) at z ~ 3, the IGM is relatively transparent. 
At shorter rest-frame wavelengths, the IGM trans¬ 
mission is significantly lower. The reason for this 
difference is that the mean-free-path of a L yC p ho¬ 
ton a t z ~ 3 is l m f p ~ 100 Mpc (IProchask a et al.1 
12009 1. As a photon travels that distance, it red- 
shifts by about 10%. Therefore, a photon with 
A res t = 829 A at 2 ~ 3 has to travel 1 l m fp to red- 
shift beyond the Lyman limit so that it is no longer 
subject to bound-free absorption from hydrogen. A 
photon of A rest = 760 A has to travel 2 l m f p (Note, 
we are ignoring the much smaller effects of decreas¬ 
ing cross section with decreasing wavelength and 
the changing mean-free-path with redshift to illus¬ 
trate the simple point.) In summary, surveys that 
probe the LyC of a large number of galax ies at 
an op timal redshift for the filter (e.g. iNestor et all 
l201lL at z ~ 3.09) will likely have a higher detec¬ 
tion rate than a survey of field galaxies at a vari¬ 
ety of redshifts. Of course, if the observations are 
done spectroscopically, then one can analyze many 
galaxies at a variety of redshifts with good sensi¬ 
tivity near the Lyman limit. 

• Filter width: Because the mean-free-path of a LyC 
photon in the IGM at z ~ 3 corresponds to ob¬ 
served 330 A (see above), the filter width can af¬ 
fect the level of signal detected. Filters signifi¬ 
cantly smaller than the mean-free-path (and near 
the Lyman limit) will probe regions of the LyC with 
high transmission. Broadband filters however, will 
probe wavelengths with little or no transmission 
at the blue end of the filter, as the typical filter 
(e.g. the Johnson U filter) is twice the width of 
the mean free path. In fact, a wide filter will of¬ 
ten be adding noise at the short wavelength end, 
but no signal. Therefore, if two images reach the 
same depth in AB magnitudes but one uses a nar¬ 
row band (near the Lyman limit) and the other a 
broad band, the narrow band will be more sensitive 
to escaping LyC. 

All of these factors must be carefully considered when 
comparing candidate detection rates and contamination 
rates. If the detection rate is not significantly higher 
than the expected contamination rate, then no detection 
should be claimed. 

INestor et all (1201 111 carefully calculate the contamina¬ 
tion rate using the surface density of objects at the same 
wavelength as the LyC filter, and using a contamination 
radius (1.2") equal to the radius used to match the target 
catalog to the LyC filter catalog. When doing so, they 
determine a contamination rate that would account for 
about half of the detected brighter Lyman break-selected 
galaxies, and ~ 16% of the fainter, Lya-selected galaxies 


(LAEs). Our results in this paper are consistent with the 
high contamination rate for the bright er LBGs. 

For comparison, llwata et~ahl (120091) covered a much 
wider area, but 0.6 magnitudes shallower in depth. Be¬ 
cause of this, they were sensitive only to the escaping 
LyC in more luminous galaxies. They detect a very low 
fraction of their candidates, 8.6% (17 of 198), and es¬ 
timate that only ~ 2.37% (~ 4.7 detections) could be 
explained by foreground contamination. However, the 
contamination rate estimate was not precise. They used 


published number counts at 3000 A, not the wavelength 
of their observations (~ 3590 A) . As number counts rise 
quickly from the near-UV to optical, the counts turn out 
to be significantly different. Also, the counts they use are 
incomplete, and are therefore too low. If instead we use 
the completeness-c orrected cumulative number counts 
at ~ 3 750 A from iNonino et all (I2009D : lYanzella et al.l 
(|2010bl) . then we get a much larger foreground source 
density (~ 3.5 x 10 5 deg -2 instead of 1.0 x 10 5 deg -2 ). 
This estimate increases the expected contamination by a 
factor of 3.5 and the expected number of contaminants 
is ~ 16.6, consistent with the 17 detected sources. It is 
true that, because the number counts stated above are 
from slightly redder wavelengths, the foreground source 
counts are slightly higher than at ~ 3590 A. However, 
the difference is likely small because of the small wave¬ 
length change, and even if the contamination rate were 
increased by only a factor of 2.5 (to 5.9%, or 11.7 con¬ 
taminants), instead of 3.5, from their original estimate, 
the significance of 17 detections would only be a 1.4 a 
deviation from the expectation value. Of course, there 
could be real LyC sources in the sample, but it cannot be 
claimed with high confidence that the sample is anything 
but foregr ound sources. 

IBoutsia et al.l (I2011D were careful to only analyze 
galaxies at redshifts where the Lyman limit was close to 
the red edge of transmission in the filter (3.27 < z < 3.35, 
when the transmission edge corresponds to z ~ 3.27). 
However, their sample was limited to relatively luminous 
galaxies (10 of 11 galaxies have L > L*), where the aver¬ 
age escape fraction appears to be quite low. In addition, 
the sample is small, and it is very possible, given the 
rare LyC detections of luminous galaxies, that had they 
sampled significantly larger samples, they would have de¬ 
tected some galaxies. Nestor et al. only detected 20% 
of luminous LBGs with deeper data (ler = 29.1 AB in 
the narrow-band filter, which is more sensitive to LyC 
vs. 2 9.3 AB in the broadband filter), so a null result for 
IBoutsia et al.l (1201 ll) is n ot unexpected. 

lYanzella et al.1 (l2010al) had large samples and were 
careful about contamination. However, the sample is 
very bright (so a low percentage of detections is already 
expected) and many of the galaxies are at much higher 
redshift than the most sensitive redshifts. Specifically, 
less than half of the non-AGN (58 of 128) are at z < 3.6, 
such that the longest wavelength of transmission in the 
filter (~ 3900 A) is sensitive to the galaxies’ LyC emission 
withi n one mean-free-path of the ( l m fp{z = 3.6) ~ 50 
Mpc, IProchaska et alJl2009f ) . Thus, most of these galax¬ 
ies would not be detectable with significant escape frac¬ 
tions. In addition, the escape fraction limits for the 
galaxies at z ~ 3.5 — 3.6 are significantly less sensitive 
than those at z ~ 3.4, the lowest redshift probed by 
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their filter . The refore , it is perhaps not surprising that 
lYanzella et al.l (l2010af) do not have a significant number 
of detections (after using high resolution Hubble resolu¬ 
tion to minimize foreground contamination). 

In summary, all of these surveys at z ~ 3 are consis¬ 
tent with the interpretation that high escape fractions 
in bright, L > 0.5 L*, galaxies are rare. Our results in 
thi s paper ar e cons istent with this scenario. In addition, 
INestor et al.l (12013H suggest it is more common to have 
high escape fractions from faint Lya emitters. In a sub¬ 
sequent paper, we will investigate escape fractions and 
foreground contamination in the fainter sample. 


6. CONCLUSIONS 

In total, there were six bright, 2 ~ 3.1, LBGs with 
candidate LyC emission i dentified via detect ion in the 
Keck/LRIS NB3640 filter (|Nestor et al.l 120 fTI) . Here we 
have presented follow-up HST and Keck/NIRSPEC ob¬ 
servations of five of the six candidates. In each case, the 
candidate LyC emission is associated with one of multi¬ 
ple clumps or, in the case of D17, offset from the primary 
galaxy but not obviously associated with another galaxy. 
Therefore, spectroscopic follow-up was required to con¬ 
firm the redshifts of the clumps that were emitting the 
candidate LyC emission. 

One of the five target galaxies (aug96ml6) was found 
to have been misidentified as a galaxy at z > 3.06 from 
the Keck/LRIS optical spectrum. The actual redshift of 
the galaxy detected in the NB3640 filter is unknown so 
it should not have been included in the original sample 
as a LyC-emitting candidate. Therefore, there are four 
galaxies with z > 3.06 with candidate LyC emission. 
The Keck/NIRSPEC spectroscopy revealed that two of 
the remaining four candidates were found to be con¬ 
taminated by galaxies at lower redshifts (C49_NW and 
MD32_NE). A redshift for the component to the NW o f 
MD46, MD46_NW, could not be confirmed. INestor et all 
(l20ll suggest that extended Lya emission in the di¬ 
rection of MD46_NW indicates that the galaxy is likely 
at the same redshift as the brighter counterpart, MD46. 
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However, the NIRSPEC limit on [O ill]A5007 is quite 
low, suggesting that it must have a significantly lower 
L[oiii]/L uv ratio than MD46 if MD46_NW is indeed 
at the same redshift as MD46. This remains the most 
promising LyC-emitting candidate of the bright galaxies 
in this paper. Finally, the purported LyC emission from 
D17 is offset by ~ 1". There is not a counterpart galaxy 
at the location of the emission down to F814W ~ 27.3 
(AB). The low S/N emission seen from the ground is not 
seen in the F336W image, however, varying IGM opaci¬ 
ties may account for this difference. 

In summary, we find that the contamination rate of 
this bright, L > 0.5L* sample of galaxies is at least 
50% (2 of 4) and we do not confirm any definitive LyC 
emitters in the sample. This is consistent with the con ¬ 
tamination r a tes es timated in INestor et al.l (1201 11120131) . 
INestor et~ a l. (2013) find significantly lower contamina¬ 
tion rates in their faint, Lya-selected sample, suggest¬ 
ing that these types of galaxies are have far higher av¬ 
erage escape fractions. In Section [5] we show that the 
other z ~ 3 surveys for LyC are generally consistent with 
this picture, once varying depth, sample luminosity, and 
contamination rate are taken into account. Because the 
fainter, Lya-selected samples are expected to have higher 
escape fractions and lower contamination rates, in a sub¬ 
sequent paper, we will present results of additional HST 
imaging and Keck spectroscopy of this promising popu¬ 
lation. 
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